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Abstract. High resolution spectra data of red clump stars towards the NGP have been obtained with the high 
resolution spectrograph Elodie at OHP for Tycho-2 selected stars. Combined with Hipparcos local analogues, we 
determine both the gravitational force law perpendicular to the Galactic plane, and the total surface mass density 
and thickness of the Galactic disk. The surface mass density of the Galactic disk within 800 pc derived from this 
analysis is E {\z\ <800 pc) = 76 Mqpc”^ and, removing the dark halo contribution, the total disk mass density 
is Eq = 67 M 0 pc“^ at solar radius. The thickness of the total disk mass distribution is dynamically measured 
for the first time and is found to be 3901^20 relative agreement with the old stellar disk scaleheight. All 

dynamical evidences concerning the structure of the disk (its local volume density -i.e. the Oort limit-, its surface 
density and its thickness) are compatible with our current knowledge of the corresponding stellar disk properties. 
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1. Introduction 

We present a new dynamical determination of the Galactic 
potential and Kz force perpendicular to the Galactic 
plane. We measure the surface mass density of the 
Galactic disk (i.e. the total amount of disk mass in a 
column perpendicular to the Galactic plane) and also 
the thickness of the vertical mass distribution of the 
disk at the Sun. They are fundamental parameters of 
the Galaxy, determining the disk contribution to the 
Galactic potential and to the rotation curve. They also 
give valuable constraints to the disk properties. The 
problem” is traditional in Galactic dynamics and subse¬ 
quent 
limit) 


tleiivatiun of the total mass density pp (the Quit 

las a lung histuiy beginiiig with KapLeyT~fl922)' 


an d [Port (1932). Comprehensive review may be found 
in Kerr & Lynden-Bell ( 1986|) and in Kuijken ( 19951) . 
Most recent works may be found in Creze et al. (1998a) 
and in Holmberg fc Flyiml (2000) and references within. 
They favor a low Oort limit p(z=0) = 0.076-0 .IOMqpc"^. 
Previous surface mass density determinations (including 
only disk contribution and removing the local halo con¬ 
tribution) are So= 52±13 Mqpc”^ ( Flynn fc Fuchs] 1 
Eq = 48±9M 0 pc“^, (Kuijken and Gilmore 1991). 


Flere, to investigate the vertical structure and poten¬ 
tial of the Galactic disk, we have observed a sample of red 
clump giants towards the North Galactic Pole with a high 
resolution spectrograph (Elodie at OHP). Observations 
and pr operties of the sampl e are given in the paper I of this 
series ( Soubiran et al 2002| ) where we detail the measure¬ 
ments of abundances, absolute magnitudes and radial ve¬ 
locities, and describe the sample properties. Atmospheric 
parameters and absolute magnitudes of stars have been 
obtained with the TGMET method (Katz et al. 1998) by 
comparing the observed spectra to the TGMET library 
of similar spectra of Hipparcos stars with known charac¬ 


teristics. The method is applied in Paper I (Soubiran et 
al 2002| ). The library is built from the ELODIE database 


Send offprint requests to: siebert@astro.u-strasbg.fr 
* Based on observations taken at the Observatoire de Haute 
Provence (OHP) (France), operated by the French GNRS. 


described in Prugniel & Soubiran (2001). The main ad¬ 
vantage of selecting red clump stars is that their luminos¬ 
ity function is sharply defined in a small magnitude in¬ 
terval minimizing completeness bias (Malmquist or Lutz- 
Kelker), making more accurate their corrections. This re¬ 
mote sample has been analysed combined to a local sample 
of similar Hipparcos stars in a sphere of 125 pc around the 
Sun. 

In Sect. 2, we describe the red clump sample definition 
extracted from the general sample of bright red stars ob¬ 
served at OHP (see Paper I) and the correction for the 
Lutz-Kelker bias for the associated Hipparcos subsample. 
The method to determine the potential $(z) and vertical 
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force Kz is classical but calls for some algebra because the 
NGP samples cover wide solid angles on the sky (Sect. 3). 
Results are presented in Sect. 4: we confirm with this in¬ 
dependent stellar sample previous hndings concerning the 
disk surface mass density and we present the first obser¬ 
vational constrain on the total thickness of the mass dis¬ 
tribution of the Galactic disk. 


2. The giant clump samples 


Our sample c onsists of st ars selected out of the Tycho- 
2 catalogue (ESA 


_ ^ 1997 ) and is described in Paper I 

(^oubiran et al 2002). A selection in B — V colour has been 


applied to increase the number of red clump stars with re¬ 
spect to dwarfs, subgiants, RGB and AGB stars. We have 
chosen the colour interval 0.9< B—V <1.1 that optimizes 
the detection of red clump stars. The lower B — V limit 
at 0.9 rejects the brightest main sequence stars (Mv < 5), 
some subgiants and also most giants of the blue horizontal 
branch (the ’’clump stars” with the lowest metallicities). 
The upper B — V limit at 1.1 rejects most of RGB and 
AGB stars. 

The selection criteria maximise the number of clump stars 
with respect to dwarfs and red giants along the ascend¬ 
ing branch but remain a mixture of those three types of 
objects. Our stellar parameters determination (absolute 
magnitude) allows to separate red clump stars from dwarfs 
and from other giant stars. We combine our remote sample 
with a local Hipparcos sample built with the same criteria 
based on colours and abs olut e magnitudes. This procedure 
is sketched in subsection 2.1. The distance limited sample 
is built, using the Hipparcos catalog. This sample is af¬ 
fected by the Lutz-Kelker bias, its compos ition as well as 
the bias correction are given in subsection 2.2. 


2.1. The NGP cone samples 

In this study we are interested in the bright stellar pop¬ 
ulation defined by the red clump stars allowing to probe 
stellar populations far from the mid-plane. 

Red clump stars were selected on an absolute magni¬ 
tude and a colour criteria. Those criteria are defined to 
maximize the number of stars but within a small absolute 
magnitude interval. The reason is to minimize a priori 
the Malmquist and/or Lutz-Kelker bias. Moreover, isolat¬ 
ing clump stars allows to define a well defined sample in 
term of stellar evolution. On Fig. [|, red clump stars ap¬ 
pear as an overdensity mainly in the [My, logg] and [logg, 
Teff] planes enabling us to perform a reliable selection. 
An absolute magnitude interval is set, using the selection 
of the overdensity, to retrieve a maximum of red clump 
stars and a minimum of background stars. The horizontal 
lines on the two right panels correspond to My = 0.4 and 
My =1.15 which is the best compromise between the total 
number of red clump stars and the number of background 
objects in the resulting sample. Gircles on the four panels 
show stars that do not belong to the overdensity (i.e. the 
red clump), but belong to the color and absolute magni¬ 


tude selection. The number of stars in the sample rises 
from 221 (overdensity selection) to 232 (absolute magni¬ 
tude criteria) showing that a small fraction (< 5%) of 
background stars are included. 

The histogram on Fig. ^ draws the presence of the red 
clump as a narrow peak at Mv=0.8. The two vertical lines 
at Mv = 0.4 and 1.15 are our selection interval. Only a 
few other type of giants contribute to the background. 

This sample (two fields towards the NGP) is separated 
in four subsamples, two distance-complete samples and 
two magnitude complete samples, one for each field. The 
resulting number of stars in each subsample is given in 
Tab. 1^, to be compared to the 537 similar stars extracted 
from the Hipparcos catalogue (next section). 



Fig. 2. Absolute magnitude histogram of the observed 
NGP samples. The vertical lines define the selection of 
the red clump stars around the central peak. 


2.2. The Hipparcos sphere sample 

A local sample of clu mp giants h as been drawn from the 
Hipparcos catalogue ( ESA 1997| ) which provides B and 
V magnitudes, parallaxes and proper motions. This sam¬ 


ple is corrected from th e Lutz-Kelker bias using Turon 
[Lacarrieu fc Crez6 (1977) method which is outlined in sec¬ 


tion 2.2.1 . The selection criteria used for extracting red 
clump giants is transposed to the Hipparcos sample in or¬ 
der to ensure the homogeneity of the different samples. 
The sample was completed with radial velocities from the 
literature, the full de scription of its construction is out¬ 
lined in section 2.2.3. 


2.2.1. Lutz-Kelker bias and corrections 

Gonversion of parallaxes into distances introduces a bias, 
stars appearing statistically further away than they really 


are. Since the proposed correction of this bias by Lutz & 


Kelkei 

(1973 

), other methods have been proposed ( 

Turon- 

Lacarrieu & Greze 

1977; 

Arenou & Lur^ 

1999|). To adjust 
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Fig. 1. Properties of the observed stars. Horizontal dashed lines indicate our selection procedure for the red giant 
clump sample. Circles denote stars belonging to our My selection that are probably not clump stars (see text). 

Table 1. Description of the north Galactic pole samples 


Field 

dlow 

dup 

number of stars 

surface (square deg) 

Field 1 

226 

616 

85 

309.4 

Field 2 

226 

405 

49 

410.1 

Field 

V low 

V up 

number of stars 

surface (square deg) 

Field 1 

7.2 

10.1 

128 

309.4 

Filed 2 

7.2 

9.2 

73 

410.1 


the overall distribution and also the individual values, we where AM is the corrective term. It follows that 


use the process of Puron-Lacarrieu & Creze (1977) that 
allows us to correct the magnitude bias according to the 
precise s hape of the LF. Th e same method has been ap¬ 
plied by Girardi et al. ( 1998 ) to analyse Hipparcos giants. 
Here we derive the formal correction for distances. The 
later corrections enable us to defined a volume limited 
sample, free from Lutz-Kelker bias. 

Relating the observed absolute magnitude Mobs to the 
true absolute magnitude Mtrue, we can write the unbiased 
estimator of Mtrue as: 

< Mobs + AM > = Mtrue) 


<AM> = < Mtrue- Mobs > 

^ /(Mtrue - Mobs)diP(MtruelMobs) 

f diP(MtruelMobs) 


( 1 ) 


where "P (Mtrue I Mobs) is the conditional probability that 
the absolute magnitude is Mtrue if the observed one is 
Mobs- We derive dP(Mtrue|Mobs) which is proportional 
to P(7rtrueKobs)d7rtrue ( |Turon-Lacarrieu fc Gre^ 1977). 
P(7J‘true|7robs) is the conditional probability that TTobs being 
the observed value, the real parallax is TTtrue- Under the 
assumption of gaussian errors on the parallaxes and for a 
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uniform density of stars, this term is given by 


rue I'^obs) 


'^true^('^true) 

!— 0 XP 

V ^TTCTjrtrue 


('^obs ^^true) 


2(7 2 

”true 


where $(Mtrue) is the luminosity function of the studied 
tracer population. We assume a gaussian distribution for 
red clump stars with mean <M> and standard deviation 
anc- Therefore ^{M) writes 


$(M) = oq + aiM + a 2 M^ 


+ 


A 


V^CTRC 


exp 


(M- <M >)2 




( 2 ) 


where the second order polynomial represents the contri¬ 
bution of other giants. 

Calling 


^TCobs (’^true) — 


^true^(-^true)_ 

- 6Xp 

V ^'^^TTtrue 


("^obs ’^true) 


2(7 


TT 


2 

true 


Mtrue « 5 l0g(7rtrue) + TO - 5 

TO being the apparent magnitude of the star (noting that 
we are dominated by errors on the parallaxes and that 
errors on apparent magnitude are negligible), and 


(Mt,,e-Mobs) =51ogf^') 

\ '^obs / 

the correction given by Eq. ^ rewrites 

/o°° G'^obs (7rtrue)5 log dTTt, 


<AM> = 


/o (^iTobs ('^true ) dTTtrue 


(3) 


Using the same approach, we derive the mean correction 
to be applied to the distances. This correction is given by 


<Ad> = 






5 ('^true) dTTtrue 


(4) 


Equations and ^ allows us to correct on a statistical 
basis the individual values of the distance and absolute 
magnitude, given the luminosity function of our sample. 


2.2.2. Red Clump luminosity function 

To correct the observed magnitudes from Lutz-Kelker 
bias, we must determine the luminosity function of the 
tracer population. This luminosity function is obtained 
using a nearby sample built out of the Hipparcos cata¬ 
logue. We select all stars closer than 50 pc (i.e. parallaxes 
> 20 mas), with absolute magnitudes in the range [0.2,1.3] 
and with B — V colours in the range [0.9,1.1]. This selec¬ 
tion gives a distribution that can be considered as unbi¬ 
ased from errors on parallaxes. The observed distribution 
of absolute magnitudes is therefore close to the true dis¬ 
tribution and we estimate it by a least-square fit of the 
function given by Eq. |^. 


Figure || shows the result of the least-square fit and 
the observed distribution. The local clump parameters are 
< Mv >= 0.74, (TMv = 0-25 and A = 13.61. The other 
giants, modelled by the second order polynomial ( best fit 
parameters uq ~ 10“"^, ai ^ 10“® and 02 = 0.0531) give a 
contribution to the sample smaller than 2%. 



Fig. 3. Distribution of absolute magnitudes for the clump 
stars closer than 50 pc from the sun. The observed distri¬ 
bution has been fitted using Eq. || giving < > = 0.74 

and ctMv = 0-25 for the local clump. The dashed line 
shows the contribution of field giants to the sample. 


2.2.3. Construction of the sample 


The local sphere sample is built using Hipparcos data 
where we select all stars with parallaxes larger than 5 mas 
(i.e. stars closer than 200 pc) in the same colour interval 
as the two NGP samples. 

Then we correct the sample from Lutz-Kelker bias and 
we compare both corrections, first using the absolute mag¬ 
nitude correction (Eq. ^) and second using the distance 
correction (Eq. The mean absolute magnitude correc¬ 
tion of the sample is —0.09 mag with a standard deviation 
equal to 0.13 mag, and for the distances the mean correc¬ 
tion is —5.7pc with a standard deviation equal to 23.5 pc. 
Figure]^ shows the effect of the Lutz-Kelker bias on the lu¬ 
minosity function before and after correction. The contin¬ 
uous line gives the observed distribution in absolute mag¬ 
nitude, whereas the dashed line shows the corrected lumi¬ 
nosity function. The sample defined by stars closer than 
125 pc (537 stars) was then established using the corrected 
distances. The consistency has been checked by comparing 
the distance obtained after correction with Eq. to the 
distance derived using the absolute magnitude correction. 
The difference in the two distances being of the order of 
0.2 pc for the sample, this ensures us of the reliability of 
the correction. 

We use the Barbier-Brossat &: Figon ( 2000| ) cata¬ 
logue that provides us with the radial velocities. We com- 
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Table 2. Distribution of the quality index for radial ve¬ 


locity of the stars of our sample in Barbier-Brossat & 
Figoi£ 200C). 


Quality Index 


GVr 

(km/s) 

number of stars 

A 


GVr 

< 2.5 

74 

B 

2.5 < 

GVr 

< 5.0 

170 

C 

5.0 < 

GVr 

< 10.0 

23 

D 

10.0 < 

GVr 


8 


plete the sample with radial velocities from the Simbad 
database when measurements are existing. This leads to 
a nearly complete (98%) sample of 526 stars (out of 
537 stars). Among those 526 stars, 27 have radial veloci¬ 
ties from Simbad whereas 499 have radial velocities from 
Barbier-Brossat fc Figon ( ^000 ). 

Measurements errors for the radial velocities are given 
in two distinct ways in [Barbier-Brossat fc Figon] ( 2000] ), 
either a value for avr is given when available (224 stars) 
or a quality index (A to D) for stars whose radial ve¬ 
locities are obtained from a prism objective survey (275 
stars). Figure ^ shows the distribution of errors of the ra¬ 
dial velocity for stars with known avr while Tab. ^ gives 
the distribution of the quality index for the stars with no 

Gy r. 



avr (km/s) 

Fig. 4. Distribution of standard errors for the stars of the 


Hipparcos sample with given avr in Barbier-Brossat & 
Figon|~} 2000[ ) (224 stars). 6 stars have errors larger than 5 
km/s and are not within the limits of the plot. 



Fig. 5. Observed luminosity distribution (continuous line) 
versus the Lutz-Kelker bias corrected luminosity distribu¬ 
tion (dashed line). 


assumption is that the stellar motion separates in vertical 
and horizontal motions (separable potential), so the dis¬ 
tribution function f{z,w) of a tracer population depends 
only on the potential of the disk mass and on the kinetic 
energy of stars. For instance, it may be written as: 


f{z,w) = > 


P*(0) 


N 

\/^o 

1 — 1 ^ 


exp 


4>(z) -b \ 


(5) 


where N is the number of isothermal disks, the popu¬ 
lation having the vertical velocity dispersion a^i and the 
local density pi(0). The vertical density of a stationary 
population is just: 


N N 

i=l i=l 


This second assumption is valid as long as the analysed 
sample covers a restricted range of distances close to the 
Galactic plane like the bulk of our NGP stars between 
162 pc and 870 pc. At higher distances, the potential can¬ 
not be considered as separable in r and z and the in¬ 
clination of the velocity ellipsoid must be included in the 
analysis: see for instance Kuijken and Gilmore (1989) and 
^tatlei ( 1989 ). 


The vertical potential ( Kuijken and Gilmore 1989 ) is 
written 

$(z) = 27rG(So(v^z2 + D^-D)+ p^sz^) (7) 


3. The force perpendicular to the Galactic plane 
and the disk mass distribution 

Our first assumption is that stellar population distribu¬ 
tions are in a stationary state with well mixed distribu¬ 
tion functions according to coordinate z (distance from 
the Galactic plane) and w (vertical velocity). The second 


a 3-parameter {D, Eq) Pes) representation to study the 
shape of the vertical potential. Another approach con¬ 
sists in building a stellar disc mass model and in con¬ 
sidering that supplementary mass (unidentified mass) is 
proportional to the disc mass model or to one of its sub¬ 
component (see for instance Bahcall (1984]); Holmberg & 
Flynn] ( pOOO]) ). 
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The total vertical mass distribution is deduced from 
$(z) through the Poisson equation and discussed in many 
papers (con Hoernei (1960); for extensive references see 
Creze et al. ( |1998a ); Holmberg & Flynn ( 2000D , they also 
discuss the horizontal potential contribution to the dy¬ 
namical estimate of the disk mass). 

We derive the potential parameters as well as the stel¬ 
lar sample distributions (through the pi and ) by using 
a maximum likelihood method. The maximum likelihood 
method allows the ht without binning of data. It allows a 
non-biased estimate of parameters even when the number 
of objects in bins would be small and the fluctuations are 
not gaussian and the use of a least square fitting would 
not be justihed. The oretical bases for this sch eme can be 
foun d for in stance in Kendall fc Stuart ( 1973|) or in Eadie 
et all (|1971D . 


We set f*{xi) the probability function of observables 
Xi whose detailed expression is given below. 

The logarithm of the likelihood is then defined as: 


logL = ^log/*(a;i), 


( 8 ) 


i=l 


where n* is the number of red clump stars in the whole 
sample. We have separated the fields in subsamples ac¬ 
cording to distance and/or magnitude completeness. Then 
each subsample requires a specific treatment and f* has 
a different detailed algebrical expressions. It may be how¬ 
ever written as: 


r = 


Al 1 -l- Ad2 -t- Ads 

Afi + A// -l- A/3 


where the Adi are the distribution function in the volume 
Vi expressed in the used observable variables, and the A/) 
are their normalizations over each volume V). The Adi and 
A/i are explicited in the following subsections according to 
the type of completeness. 

We have determined the model parameters using dif¬ 
ferent observables (z, r, 5, w, TOv, or My) in order to find 
the most efficient and accurate way to analyse the samples. 
Model with a simple limit (z-distance) have a simple al¬ 
gebra and are easier to develop. Models using supplemen¬ 
tary stars (r-distance or magnitude completeness) would 
be more accurate, they have also more complicated alge¬ 
bra. We compare the numerical consistency of the various 
methods and check the most efficient and accurate ones. 

Observations are available in three volumes V^, the 
two cones towards the NGP {k = 1, 2) and the local sphere 
(fc = 3). The number density of the stars is, according to 
the dynamical model: 


dN{x,y,z,w) = 


f{z, w)dxdydzAw 
0 if {x,y,z,w) is outside 


or in galactic coordinates: 

/(rsin(6), w) cos(5) drdfedZdrc 


dA^(r, /, b, w) = 


Due to the extended angular size of the observed NGP 
fields, the dependence on Galactic latitude b cannot be 
neglected and, by comparison to traditional studies, makes 
the algebra more tedious. 

3.1. Local sphere 

The sample extracted from Hipparcos data and corrected 
for Lutz-Kelker bias defines a volume that is complete in 
distance. This volume 14=3 is a sphere of radius Rs = 
125 pc. 

The stellar density according to variables z and w is 
modelled as: 

dJV(z w)= I^ 

^ \0 if |z| > Rs 

We have = dA/ 3 /(dzdw) and 

A/'s = [ dN = 'Si [ pi{z) TT {R^^ - z'^) dz 
JVi Jz 

3.2. NGP: completeness in z-distances 

We have used the largest cone samples towards the NGP 
that are complete in z height above the Galactic plane. 
For held 1, the volume Vi is limited in z between 226 and 
391 pc (309.4 square degree held) and for held 2, between 
226 and 579pc (410.1 square degree held). 

We have with fc = 1, 2: 

r / f I f L/ d 1 dl] f(z,w)z'^dzdw 

dA4(2,u;) = M ^ 

[0 outside of the volume 

where 6iow and 6up are the lower and upper latitude bound¬ 
aries at hxed longitude. We have Mk = dA4/(dzdw), and 


A4= / dN = 

JVk 


1 


n bu 


2(sin 6)2 


dM J Piiz) z^dz 


Integrals can be separated since the I and 6 contours do 
not depend on z, according to the dehnition of the sub¬ 
samples. 

3.3. NGP: completeness in r-distances 

To increase the number of stars used to constrain the 
model, we dehne complete volumes in distance from the 
sun, r. For volume Vi, r range from 226 to 405 pc and for 
volume V 2 , r range from 226 to 616pc. 

We have with fc = 1, 2: 


(Wk{b,r,w) = 


j dN = AZ(6) /(z, w) cos{b)drdbdw 


0 outside of the volume Vj, 


0 if (r, I, 6, w) is outside Vfc 


we dehne Mk = dA4/(drd6dw), and 

Afk = [ dA4 = Si / Al{b) f pi(z) cos(6)drd6 
Jvi Jb Jr 



































Siebert et al.: The surface mass density in the Galactic plane 


7 


In the previous equations, since the distribution func¬ 
tion / is independent of I, the Al{b) term, which is the 
longitude interval at latitude 6 of a circular field, has a 
closed form given by solving r • = cos(a), where a is the 

radius of the field, the unit vector pointing towards the 
center of the field and r a unit vector. With 1^ and he the 
coordinates of the center of the field and a its radius, we 
obtain the expression valid for the field 1: 


Al{b) = 2 arccos 


cos(a) — sin(6) sin(6c) 
cos(&) cos(6c) 


For the field 2, Al{h) is 27r since the field is centred to¬ 
ward the Galactic pole, unless the latitude b corresponds 
to the removed area of the Coma Berenices cluster or to 
the field 1 overlapping region. Since these removed region 
are circular, we have: 


Al{b) = 2 tt- AcomJib) - Afield il{b) ■ 


4. Results 

Modelling. The modelling of the samples consists in repro¬ 
ducing the distribution of the observed apparent magni¬ 
tudes, absolute magnitudes, distances and vertical veloc¬ 
ities by adjusting the three vertical potential parameters 
(the total surface mass density Soj the mass scale height 
D, and the halo local effective mass density peff) as well 
as the stellar density and velocity distributions (through 
the velocity dispersions and the respective contribution of 
each isothermal component). Using the four methods pre¬ 
viously described, the model parameters are adjusted by 
a maximum likelihood. The maximum likelihood method 
avoids bias that would be introduced by a binning of data 
and by the small size of our samples. As expected, the 
accuracy is better with methods using more stars and/or 
more distant stars, while the comparison of methods al¬ 
lows to check their respective robustness. Following discus¬ 
sions are based on the results given by methods number 3 
and 4. 


3.4. NGP: apparent magnitude completeness 

Since the NGP samples have been selected by apparent 
magnitude, it is more efficient, in terms of the number of 
stars used, to keep this criterion. Modelling the distribu¬ 
tion of apparent magnitudes of selected clump stars allows 
the use of 201 stars from 160 to 870 pc from the Galactic 
plane. 

With fc = 1,2, we have: 

Aik = dA4(&, rriv, w)/{dbdrrivdw) = 

( = Al{b) f{r sin(6), w)^{My) cos(6) dr 
[or =0 if ruv is outside the limits given in Table ^ 

where Al{b) has been previously defined. We have also: 


A4 = S. [Al{b) 
Jb 


Pi{r sin(5))$(Mv) r^ cos(6) drd6dmv 


3.5. NGP: apparent magnitude completeness (with 
absolute magnitudes and distances) 

We use the same selection criterion in apparent magni¬ 
tudes, with 200 stars in the range 200-750 pc from the 
Galactic plane. 

Here we model both the distributions of absolute magni¬ 
tudes and distances expecting a more constraining use of 
the observables My 
We have with A: = 1,2 

Mk = dNk[z,b,w,MY)/{drdbdwdMy) = 

J = Al(b)f(z, w)d?(Mv)r^ cos(b) 

[or =0 if rriv is outside the limits given in Table ^ 

and we have: 

Afk = '^i fAl{b) f [ pi{rsm{b))^{My)r^cos{b)drdMydb 
Jb Jr J Mv 


Peg. The parameter peff is introduced to model the lo¬ 
cal contribution of a spherical massive halo. A realistic 
range for this param eter is 0 to 0.02 Mq pc“^ ( Kuijken 
land Gilmore 198S| ). However, this halo contribution to 
the vertical potential is quadratic in x (at small z) and 
cannot be clearly distinguished from the disk contribu¬ 
tion (the remaining terms in Eq. ^ that is also quadratic 
for small z. For these reasons, peff is set at 0.01 Mq pc“^. 


Adjusted parameters. The best decomposition of the ob¬ 
served stellar distributions in kinematic components, us¬ 
ing methods 3 or 4, is obtained with three isothermal com¬ 
ponents with densities (22%, 59% and 19%) and respec¬ 
tive vertical velocity dispersions Uu; (8.5, 11.4 and 31.4) 
kms“^. A two-component decomposition gives nearly the 
same likelihood with respective densities (78% and 22%) 
and vertical velocity dispersions (12.0 and 30.4) kms“^. 
This result can be compared to the kinematic decomposi¬ 
tion obtained in paper I ( ^oubiran et al 2002| ) using the 
3D velocity components and the metallicity. 

The maximum of likelihood gives also the following 
result: with method 3 {D, Eo)= (492pc, 72 Mq pc“^) and 
with method 4 {D, Eo)= (1500pc, 193 Mqpc”^). At a 1- 
cr error level methods 3 and 4 give only a lower limit for 
D respectively 227pc and 400 pc and no upper limit. The 
So values are strongly correlated with D (see Figure 
We must consider that our sample puts lower limits to 
{D, Eg) parameters and a sample with higher z-distances 
would be needed to put upper limits. 

If D were known independently, the relative accuracy 
on So would be 12 %, while considering the seven model 
parameters (including dispersions and relative density of 
the two stellar components, D and Eq), we find a large 
correlation between the two parameters D and Eq. We 
may note that we obtain the first direct observational 
constraint on the thickness of the total disk mass. We 
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also obtain a constraint on the total local mass density 
Ptotai(z=0)=0.08±0.01 M 0 pc“^ (see below). 



Fig. 6. 1, 2 and 3 -(t error contours for Eq and D solutions 
obtained with method 4 (continuous curves). Using also 
t he recent Oort limit deter mination 0.0102 Mq pc“^ from 
( Holmberg fc Flynn |2000| ) gives better constraints (1-a 
error level: dashed curve). Dashed-dotted lines show solu¬ 
tions for adopted local mass density ptotai values between 
0.06 to 0 . 2 Mqpc“^. 


Comment about D. In some previous vertical Galactic 
force determinations, simplifying assumptions have been 
adopted concerning the exact shape of the vertical 


mass distribution. For instance, Kuijken and Gilmore 
(1989) modelled the potential using Eq. |^, but assuming 
11 = 180 pc. They argue that their result is not sensitive 
to a small change on D since the majority of their sample 
ranges beyond 500 pc from the Galactic plane. However 
since a significant fraction of their sample ranges between 
200 and 500 pc, we suspect, that with a drastic change of 
D at ^500 pc, their data would be fitted with a higher 
Eq value. Other previous determinations of the total sur¬ 
face mass density, based on one parameter models for the 
potential, could suffer a similar bias. 


Table 3. Eq and D solutions for the potential within l-cr 
error and associated total local mass density and column 
density. 


Eo 

M0pc"^ 

D 

pc 

Pea 

Mopc”^ 

ptotal (0) 

M0pc“® 

E(800pc) 

M0pc“^ 

E(l.l kpc) 

M0pc“^ 


64. 

90. 

193. 

400 

637 

1500 

0.01 

0.01 

0.01 

0.090 

0.081 

0.074 

73.2 

86.4 

106.8 

82.2 

99.9 

136.1 


Total local mass density ptotai(z=0). A set of Eq and D so¬ 
lutions within l-cr error are given in Table 2 with their cor¬ 
responding total local mass volume density ptotai(^ = 0 ) 
(according to Eq. ^ the local mass density is given by 
Ptotai(z = 0) = Eo/(2£)) -I- peff). We also give the total 
mass density within 800 pc from the Galactic plane in¬ 
cluding the halo contribution parametrized with peff (our 
surface mass density determination applies only up to 
~ 800 pc, the limit of our sample). The acceptable range 
of values for Eq and D is large, resulting in a wide range 
of possible mass distributions. This can be seen on the left 
panel of figure ^ where are plotted the solutions given 
in Tab 2 (recall that is proportional to E(z)) . 

To disentangle the different solutions, we need more 
observations in the range 750 pc to 1.1 kpc. This is il¬ 
lustrated on the right panel of Fig. by the predicted 
clump giant density (the observable quantity), where we 
can see that the solutions quoted in Tab. || differ signifi¬ 
cantly above 750 pc. 

We may remark that the total local mass density 
Ptotai(z = 0) given in Table 2 ranges below 0.10 M 0 pc“^ 
and is in agreement with recent and independent determi¬ 
nations based on Hipparcos data, compatible with local 
mass density of known matter. 


Additional constraint from ptotai(-z = 0): the Oort limit. 
In Table 2, the different solutions differ by their pre¬ 
dicted total volume mass density, from 0.07 to 0.09 M 0 
pc“^ but remain in agreement with the recent measures of 
/Ototai(-z=0) based on Hipparcos data (|Greze et al. 1998a , b|; 
Pham 1998; Holmberg & Flynn 2000| ) ranging from 0.076 
to 0.102 M 0 pc“^. 

The data set of clump giants that we have used does 
not allow to constrain more efficiently the local volume 
mass density ptotai(^ = 0). Most of the constraints would 
come from the local sample within 125 pc, but due to the 
high velocity dispersions of the red giant stars, their den¬ 
sity distribution is nearly uniform within 125 pc and the 
bending of the vertical density distribution due to the po¬ 
tential cannot be measured significantly. We notice that 
the analysis of red clump stars in a sphere of 250 pc would 
help to constrain the local volume density since they have 
distances measured by Hipparcos but they have no mea¬ 
sured or published radial velocities. This would give a new 
independent determination, however the amplitude of the 
Lutz-Kelker or Malmquist bias will be large and certainly 
difficult to model properly. 

So we consider the recent results based on Hipparcos 
data of the local volume mass density |Otot ai (-g=0) - 
Diffe rences on ptotai(-z= 0) betw een Greze et al. ( 1998aD 
and Holmberg & Flynn ( 2000|) are within l-cr and are 
explained by Holmberg & Flynn ( 2000| ) as due to dif¬ 
ferent assumptions on the exact shape of the potential. 
preze et~ (1998a) assume a quadratic potential close 
to the Galactic plane while Holmberg fc Flylii^ (2000) 
suppose that the vertical potential is proportional to a 
more realistic model of the vertical stellar disk distri- 
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Fig. 7. Left : vertical force in (km/s)^/pc for the three solutions within l-cr presented in Tab. 2. Right : associated 
stellar density in log (Mq pc“^). 


bution. Admitting this more realistic hypothesis for the 
shape of the potential, we consider that Hipparcos data 
gives ptotai(^ = 0) = 0.102 ± 0.010 M© pc“^ 

( Holmberg fc Fly^ 2000 ) and combining this result and 
our result obtained from distant red clump giants, we de¬ 
duce stronger constraints on Eg and D (Figure |[ dashed 
contour draws the l-cr limit). We deduce that the scale 
height D is relatively large, SOOl^^oP*^’ 
total surface mass density Eg = Mqpc"^, and the 
mass within 800pc is E(|z| <800pc) = 76^12 AloPC”^. 
The quantity E(|z| <1100pc) =85^13M©pc“^ can be di¬ 


rectly com pared to the result obtained by Kuijken and 
Gilmoi[e~| ( 1991 ) : E(|z| <1100pc) = 71±6M0pc“^. The 
two results are compatible in the Icr limit but our error 
bars are very large due to the poor constraint on D. The 
agreement between the maximum likelihood solution and 
the data is shown in figure The NGP sample exhibits 
an offset of the w velocity towards more negative values 
than the expected -7 km/s due to solar motion. This off¬ 
set is visible at all distances, however its increase with z 
is doubtful. We have to consider that a negative offset is 
a real feature in the data because our w velocities have 
an accuracy better than 1 km/s due to the contribution 
of the ELODIE radial velocities by more than 94% in the 
considered direction. Nevertheless, this offset is not con¬ 
stant and is therefore not accounted for in the modeling. 
Table ^ shows the mean value of w, standard deviation 
and standard error on mean W for various intervals in z. 


5. Discussion 

With samples of giants in the solar neighbourhood up to 
distances of 800 pc toward the NGP, we solve for the ver¬ 
tical Galactic potential and the total disk surface mass 
density at the solar Galactic position. 

We find a vertical potential $(z) compatible at high 
(~ 400-800pc) with previous estimates (Kuijken and 


Gilmore 1991; Flynn & Fuchs 1994). We are also able 


to measure for the first time, the thickness of the total 
disk mass distribution, and find a characteristic height 
11=390 pc with a 1 -ct range from 271 pc to 720 pc. 

This scale height D corresponds to a 350 pc scale 
height of a vertically exponential density law]^ in the range 
of z-height 200 to 500 pc. Figure Js shows also the result¬ 
ing total mass density from Eq. ^ and the stellar density 
for two exponential disks (old and thick disks respectively 
with 350 pc and 750 pc scale heights and a relative density 
of 15 per cent). This measure of the mass scale height is 
in agreement with the recent determination of the stellar 
disk scale height, h^=330 pc, by Chen et al. (2001) based 


on SDSS star counts and compatible with lower estimates 


from star counts (for instance 260 ± 50 pc by Ojha et al 
( 1996D , 2 50±60pc by [Vallenari et ah ( 2000| ) or 260±90pc 
by ^ohn ( 2002| )). Below z=600pc, the contribution of the 
stellar thick disk remains always much smaller than the old 
disk, and beyond 600 pc its contribution remains smaller 
than the dark halo (see Figure ||). 

We find for the surface mass density within 800 pc from 
the Galactic plane a value of E(800pc)= 76M0pc“^ (this 
value includes the contribution from known or “seen” mat¬ 
ter -stellar and gas- and a contribution from a round mas¬ 
sive dark halo) and we find for the total disk surface mass 
density Eg= 67M0pc“^. This disk with a round dark halo 
(with a local density ~O.OlM 0 pc“^) may produce the ob¬ 
served flat rotation curve. Modellings with similar param¬ 
eters can be found in Kuijken and Gilmore (1991) and 
|Creze et M ( 1998a ). 


^ The validity of exponential density models to represent the 
vertical structure of the disk is questioned and discussed by 
^) and he explains some of the system¬ 
atic discrepancy between authors and that at low distances 
(below 500pc), a single exponential is inadequate under any 
reasonable scenario. A practical consequence is that an expo¬ 
nential fitting to star counts is valid only in restricted range of 
distances and must not be extrapolated to z=0. 


Haywood et al. (1997< 
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Fig. 8. Predicted counts for our best fit model of the vertical potential (solid line) compared to the observed distribution 
(histogram). Top : Hipparcos sample. Middle : NGP sample of radius 15°. Bottom : NGP sample of radius 10°. Left 
: cumulative distribution of z distance (Hipparcos sample) or Vj (NGP samples). Right : distribution of vertical 
velocities. The dashed and dotted lines are the contribution of each isothermal population. 


Our 1-cr error upper limit allows a maximum 
So = 114M0pc“^ and a scale height of 720 pc. It is about 
the upper limit ^ IOOMqpc"^ putting all the Galactic 
mass in the disk (maximal disk) and still expl aining the 
observed Galactic rotation curve ( Sackett 1997 ). It would 
also imply that the dark matter is within a disk of about 


1 kpc thickness. Such a fiat dark matter has been rejected 
by some observational evidences: for instance on the basis 
of the observed coupling in the stellar 3D velocity dis¬ 
tribution in the solar neighbourhood ( |Bienayrni~| |1999D 
implying a more or less round dark halo, or also on the 
non-precessing circular orbit of the Sagittarius tail around 
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Table 4. w behaviour as a function of z. 


z interval 

number of stars 

< w > 

standard deviation 

mean error 

-125 

0 

252 

-7.8 

18.4 

1.2 

0 

125 

270 

-7.9 

17.1 

1.1 

125 

250 

30 

-13.8 

11.7 

2.1 

250 

375 

46 

-12.4 

22.5 

3.3 

375 

500 

56 

-7.7 

27.9 

3.7 

500 

625 

64 

-15.5 

26.4 

3.3 

625 

750 

52 

-9.3 

31.6 

4.4 



Fig. 9. Vertical distribution of the volume mass density in 
the Galactic disk deduced from our dynamical determina¬ 
tion (continuous line). A disk vertically exponential with a 
350pc scale height (dotted line), and complemented with 
a thick disk (scale height 750 pc, and relative density of 
15 % at z=0 (dashed line). 


our Galaxy (Ibata et al. 2001) implying also that the halo 
is most likely spherical. 

We conclude that the local volume mass density, the 
surface mass density of the disk are in agreement with 
our current knowledge of the known volume and surface 
density from gas and stellar components. Moreover, the 
thickness of the disk mass density distribution is compat¬ 
ible with the thickness of the stellar old disk. 
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